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1. INTRODUCTION 

Wind energy is a renewable source of energy which does not exude greenhouse gases into the 
atmosphere [1]. Indeed, it is notorious for its cleanliness and cost-effectiveness [2]. According to the last 
intergovernmental panel on climate change (IPCC), renewable energy production will be increased to 67% by 
2050 to preclude the proliferation of global warming [3]. Low-cost energy production is needed to fulfill this 
increase [4]. Industrialization, population explosion, and urban growth, among other social processes, have 
been surrounded by technological development [5], wherein energy usage in its different forms stands out [6]. 
Today, myriads of instruments and devices that facilitate human life are powered by electricity [7]. As a 
corollary of the burgeoning electricity demand, wind energy has been capitalized on in several wind farms' 
installations worldwide [8] thanks to the propitious prospects it carries in the field of energy production [9]. 
Generally, Wind farms are erected in places where the wind level is prolific to sustain continuous electricity 
production [10], [11]. To this end, wind farms must be placed in the most appropriate location to ensure 
lucrative wind energy consumption [12]. The placement of a wind farm is a vital determinant of its power 
production. On the other side of the coin, one of the colossal hurdles that come into being when laying down a 
blueprint of a wind park is the wake problem [13], [14]. By the same token, the existence of this effect 
encumbers the uniformity of the wind speed [15]. 

In the literature, a study conducted by Mosetti and colleagues [16] sought to determine a solution to the 
wind turbine’s location. Incidentally, a genetic algorithm (GA) was used to locate the turbine in a wind farm of 
10x10 cells. By extension, Grady et al. [17] expanded on Mosetti et al. magnum opus by carrying out the objective 
function of energy and cost of production to generate ample total power production. This operation was conducted 
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by including more turbines into the design. On the other side of the coin, Gao et al. [18] employed a related form 
of the genetic method to get the optimal layout under average conditions. To corroborate its effectiveness in 
tackling realistic conditions, this method was used in Hong Kong to a genuine offshore park. In a similar spirit, 
net present value (NPV) used by Mora et al. [19] as the objective function. In the GA, integer codification was 
applied to explore feasible efficiency-improving measures. Conversely, Emmami et al. [20] regarded the total 
energy produced as the objective function. Whereas Song et al. [21] utilized an evolutionary algorithm to optimize 
the wind farm wherein he employed complex calculations. To tackle the updated bi-criteria optimization model 
with electricity production and constraint violation as the objective functions, Kusiak et al. [21] developed a 
multiobjective evolutionary strategy method. The procedure converts each constraint into a different goal 
function, the second of which is lowered to zero. Marmidis et al. [22] utilized the Montecarlo method intending 
to optimize the wind farm location. In doing so, in steady direction, steady speed wind situation, they gathered 
more accurate findings than GA. Instead, Salcedo Sanz and colleagues [23] proposed the coral reefs optimization 
(CRO) method to tackle offshore park’s construction. In the same line of reasoning, A levelized cost of energy 
(LCOE) was performed in Chen et al. study [24] as like an objective function, and the greedy method was 
employed to determine where turbines ought to be located on the farm. In Wang et al. study [25], and to improve 
power production, the wind turbine location was addressed using a differential evolution technique and a brand- 
new encoding technique. Each wind turbine was treated as a unique in the new encoding technique and the 
population of wind turbines as a whole. In addition, the study by Yang et al. [26] who worked on a Boolean code- 
based modified genetic algorithm. Lastly, Zergane’s piece of research [27] is encapsulated in the incorporation of 
the method of pseudo-random number generation. 

The primordial purpose of this study lies in creating a suitable optimization model necessary for the 
installation of wind turbines inside a square farm. To realize this objective, MATLAB was deployed to create 
and solve the optimization model, which was in turn developed with C language. Thanks to the particle swarm 
optimization (PSO) model, the solution was successfully obtained by virtue of its capacity to burgeon the 
highest power output. In light of the issue this research endeavors to demystify, this piece of research is 
tremendously significant for several reasons as it takes into consideration the variability of both directions as 
well as speed using the metaheuristic PSO. Thus, the present study provides avant-garde contributions to the 
literature wherein two novel cases: variable speed and wind direction were put under the microscope using the 
PSO method which proved to be highly effective. With that in mind, the results were compared to other 
metaheuristic methods such as GA [17], Pseudo-number generation method [27], and adapted evolutionary 
algorithm using the Boolean coding technique [26]. Hence, this study is structured in the following manner. 
The first section unveils a general explanation of the subject. It also exhibits a comprehensive review of the 
literature as well as the previously used methods pertaining to this field. The second section provides lucid 
scrutiny of the wake effect, encompasses a rigorous description of the PSO method in addition to other 
alternative methods, and presents the field data along with a brief formulation of the problem. Lastly, the results 
of the MATLAB model were compared and were thereby evaluated. 


2. METHOD 
2.1. Mathematical modelling of the wake effect 

To produce electricity, each turbine pulls out from the wind kinetic energy. When the wake effect 
occurs [28], the wind speed decreases each time the wind passes through a turbine. It affects the environment 
in two ways: It results in a reduction in wind speed [14], which ultimately dwindles the turbine’s output power 
and leads to a decrease in life span due to the increase in turbulence intensity. There exists a scalar that identifies 
the wake and indicates how the expansion of the wake occurs at a specific distance. The scalar is denoted k. In 
(1) expresses the relationship between this scalar and the parameters of the wind turbine, the height of the hub 
is denoted Z, and the coefficient of roughness is indicated by Zo, 


Z 1 
) () 


To make the theory of the wake phenomena easier to understand, the example of two turbines is used 
in Figure 1, one turbine producing the effect and the other one being the victim of the effect. In reality, one 
wind turbine can be the victim of several wake effects, which are generated by several downstream turbines. 

The total velocity deficit according to the Jensen model at the position 1 is, 


1 
Uitotar = (Ljesci wij”)? 2) 
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Figure 1. Multiple wake effect 


We include all turbines that cause the wake effect and that affect a wind turbine at a position denote i, in a set 
S (i). This set is defined, 
Vea yy 
SW = {9 py — x;| < .| (3) 


The following is the new radius taken into account when calculating the speed of a turbine influenced by the 
wake effect (where turbine j€S (i)), 


4 = kyy +7, (4) 
The expression of wind speed at turbine in position 1 will be expressed like, 
Uj = Up (1 — Uijtotar) (5) 


Up: factor of axial induction. 
To compute the velocity loss when there are several wake effects, two other factors will be taken into 
consideration: 
— A; : rotor area of turbine i 
A;;: Intersection between the area created by the rotor of turbine i and the area of wake created by 
turbine j. 


ie att (6) 


a: axial induction factor. 
Betz relations show the link between Cr (The thrust coefficient) and the axial induction factor a, 


Cr = 4a(1 — a) (7) 


The Figure 2 illustrates the intersection between the area created by the rotor of turbine i and the area of wake 

created by turbine j that generates the phenomenon. 
The expressions of the angles shown in Figure 2 are, 
Tre +5 72-7 j7 
d= 2 cos" * (Hes it ) 
2TyXij 
r j2+x; ae (8) 
@, = 2cos"+ (S—) 
27 jXij 


The Figure 2 shows that we can distinguish between three possible cases concerning the value of A;;, 
0, for xj; > 1 +7; 
2 
Aij = TY", f OY Xiz <1] — T (9) 
1 ; ; 
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From a wind turbine the power that can be pulled out must be expressed [3], 
P == CppAU? (10) 


Where: Cp the power coefficient, p: the density of air (1.225 kg/m3), U: the wind speed and A: the rotor swept 
area. 


Partial wake area 


Wake effect area from 


upstream turbine : oes 
Wind turbine (i) rotor 


area (A;) 


Wind turbine j > 


Figure 2. Intersection area [17] 


As it is proved by (10), a wind turbine i generates a power at its output which is proportional to the 
wind speed's the cube. This output power is equal to, 


0,u; < 2.3 m/s 
0.3u;3,2.3 m/s < u; < 12.8 m/s 
630,12.8m/s <u; S$ 18m/s 

0,18 m/s < u; 


P(uj) = (11) 


2.2. Metaheuristic methods 
2.2.1. The metaheuristic of genetic algorithm GA 

Genetic algorithms describe a technique that achieve the minimum and maximum of a function using 
ideas influenced by natural genetics rather than natural transformation processes [17]. Their guiding idea is the 
preservation of a sequence of population which changes with time, trying to lead population its best.The input 
sequence can be determined by chance or by heuristics. 


2.2.2. Pseudo-random number generation 
It's an iterative approach based on deterministic functions, that are frequently used to generate pseudo- 
random numbers based on recurrent patterns [27] like, 


X(N + 1) = (aX(N) + b)mod[NC + 1] (12) 


Where NC is the total amount of turbines which does not surpass NC and a ,b are positive integers. This 
indicates that the components of this series are the remainders of dividing mostly on divisor NC+1, while 
simultaneously picking non-zero integers a and b; 


2.2.3. GA involving Boolean code modification 

The assignment probability in the binary-coded GA is a constant P;= P. As a consequence, using a 
Boolean code, the order and number of turbines will be distributed evenly that optimizes the starting 
population's variety [26]. The mutation involves converting the output of even a bit to inverse. For each bit, 
the mutation site is picked at random with equal probability (p). Mutation is vital in GAs because it introduces 
new genes into the population, preventing the algorithm from reaching a premature local optimum solution. 
The classic mutation procedure, on the other hand, still attempts to balance the amount of both Is and Os there 
in code. To eliminate this detrimental consequence, different mutation rates are utilized for codes, as follows: 
The real code for such an N-bit code has c number of 1s. Whereas if the rate of mutation is p, then p 0 in the 
Boolean code represents the mutation rate at 0 and pl represents the mutation rate for 1. The following 
describes the link between p, p 0, and p 1:mutation rate for 1. T. 


_ (N-c)potcp1 
p = Sropoters (13) 
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2.2.4. Partical swarm optimization 

The metaheuristic called PSO is an optimization tool depending on individual participation. The 
concept of this optimization algorithm PSO is that we can reach a global organization of unintelligent 
individuals [29]. All particles of the swarm gradually converge to a global minimum.The principal gain of this 
metaheuristic is that it is easy to implement because it does not require many parameters. In a PSO, each 
particle in the swarm remembers its best position as well as the best position of its neighbor because each 
particle has a large memory capacity if compared with other metaheuristic methods such as the genetic 
algorithm (GA). Initially, each particle takes its place in the search space whether randomly or not. Following 
each iteration, a particle will be moved according to three components: the current velocity V;,,the greatest 
solution P, and the greatest solution collected in its vicinity P,;. The relationship between the different 
parameters of the metaheuristic (PSO) are expressed by the two (14), 


Ve = CV + C2 (pit = xt) ae €3(Pgt = x.) (14) 
Xep1 = Xp + Veg 


V, : the velocity gathered at time t, x;: the particle’s placement at time t, p; 4: the greatest previous placement 
at time t, Pg: the greatest vicinity’s previous best at time t.cy,cz and c3 are the coefficients of social/cognitive 
confidence. 


2.3. Problem formulation 

The primary purpose of this research is to increase farm productivity by selecting the ideal position 
for each turbine, a location which the wake effect is minimal. Mosetti et al. jumped to a conclusion concerning 
the cost of wind turbines as a whole [16]. This assumption shows that a single wind turbine has an unsized cost 
per year that is one, with a maximum cost decrease of one-third for each additional rotor. Therefore, the 
normalized cost is expressed as a function of turbine’s number, 


2 2, 1 -,-0.00174N? 
Cost =N(2+=e ) (15) 
N denotes the farm's turbine number. The production that a wind turbine configuration can generate is, 
Pr = Lies Lies ficPi (Ui) (16) 


Table 1 shows the incidence rate to every angle and for each wind speed, with the sum of occurrences equal to 


unity 37° fi = 1 (Where f; is the wind frequency distribution). 


Table 1. Non-uniform wind speed, changeable direction probability density wind distribution 


Up = 8 m/s Up = 12 m/s Up = 17 m/s 
k 0-360 [0-270, 280, 290, 300, 310, 320, 330, 340, 350, 360] [0-270, 280, 290, 300, 310, 320, 330, 340, 350, 360] 
f 0.0042 [0.0084, 0.0107, 0.0126, 0.0149, 0.0149, 0.0195, [0.0112, 0.0135, 0.0163, 0.0191, 0.0302, 0.0358, 
7 : 0.0149, 0.0149, 0.0126, 0.0102] 0.0307, 0.0191, 0.0163, 0.0135] 


Since it is an optimization method, it is clear that we must define an objective function. The purpose 
of this study is to achieve the greatest quantity of electrical energy at the lowest feasible cost. So, the objective 
function is expressed, 


er C 
Objective = = 


(17) 


Ts 


To give a logical comparison of data obtained from previous research on the difficulties of avoiding the wake 
effect by improving wind turbine siting and the current study based on particle swarm metaheuristics (PSO), 
we use a model similar to the one of the previous studies: it is a representation of a wind farm by 10x10 square 
cells which depict a wind farm site of 2000x2000 m as seen in Figure 3. Each wind farm cell can be configured 
with or without wind turbines. Based on this model, a single wind turbine can have 100 locations. In Figure 3, 
we can pinpoint potential wind turbine locations. Wind speed and direction are the two wind parameters. 
Formally or informally, the wind direction can be distributed. As shown in Figure 3 the wind distribution is 
done on 36 directions. The initial direction will be equal to 0°, we denote it 0). Then, to identify the abscissa 
and ordinate of the position of wind turbine i, we use the formula, 
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{ x(6;) = x(8)) cos 0; + y(@p) sin 6; (18) 


y(0;) = —x(Oo) sin 8; + y(9) cos 6; 


WNLJUJUIUUL 


Figure 3. Divided wind farm site which is considered for case study 


3. RESULTS AND DISCUSSION 

A rectangular shape ground of 2 kilometers by 2 kilometers on which turbines will be placed has been 
investigated in myriads of previous studies. In this study, we have adopted the same condition and turbine type 
in order to achieve a plausible comparison. Figure 4 depicts the suggested algorithm's flow chart. 


Input system parameters (square Wind farm 
& Wind turbine position) 
Calculate objective function (total maximum 

ower) 


Set iteration It=0 


Initialise PSO functions & parameters 
(cl=2.5; c2=0.5; wdamp=0.99; w=1 ) 


Compare objective functions (previous & 
updated) 


Is objective 
function the best 
(min)? 


Yes 


Output the best wind farm layout 


Figure 4. The suggested algorithm's flowchart 


The study offered three instances categorized according to the number of turbines, 
— Case 1: the wind blows in one direction (6; = 0 V i) and at a specific speed (for example Up = 12 m/s = 
constant). Whenever it relates to air flow, we use the direction of the ordinate vector. 
— Case 2: In this case, we conserve a constant value of speed (for example ug = 12 m/s) and we adopt the 
scenario of multiple direction. 
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— Case 3: Such as previous work, this study will employ three alternative wind speed values: 17 m/s, 
12 m/s, and 8 m/s. For every value of speed, the direction angle 0; can take 360 values (from 0° to 360°). 
In the three cases, the number of populations will be taken such as previous studies, which mean 30 

and 39 turbines. 


3.1. Case 1: uniform wind speed and uniform direction 

For the case of 30 turbines, the use of PSO gave a minimal objective compared to the Pseudo-random 
Number Generation as shows in Table 2. The turbines configuration is presented in Figure 5 (a) and (b). 
Comparing the results obtained for each method in case | for a population of 39 mentioned in Figure 6 (a) 
and (b), in the PSO approach Figure 6 (b), the objective function has the lowest value as the first case. The 
resulting energy is 15058 kW. 


(a) (b) 


Figure 5. Optimal turbine’s configuration by Zergane et al. (a) [27] current study and (b) using 30 turbines 


(a) (b) 


Figure 6. Optimal turbine’s configuration by Yang ef al. (a) [26] current study and (b) using 39 turbines 


Table 2. Comparison results between PSO, pseudo-random number generation and GA Boolean 
Number of turbines — Objective (10-3) — P (kW) 


Modified genetic algorithm based on Boolean code [26] 39 1.515 17,768 
current study 39 1.389 18,861 

Pseudo-random Number Generation [27] 30 1.543 14,310 
current study 30 1.463 15,058 


3.2. Case 2: Changeable direction and constant speed 

The wind direction will be incremented by 10° from 0° to 360° and for a constant wind speed, which 
is equal to 12 m/s. The algorithm proposed by this study has given significant results and the objective function 
is minimal when compared with that found in the case of the genetic algorithm. A number of turbines are 
located around the periphery of the wind farm, and the wind dispersion is uniform throughout the site, therefore, 
these turbines are directly confronted with a freewind speed, so they can provide a significant power output and 
higher than others. Figure 7 (a) and (b) depicts the optimal wind farm arrangements. The turbines of the farm 
are symmetrically positioned Figure 7 (b). To prevent the creation of the wake effect, each wind turbine is kept 
above a certain spacing from the others. Since the wake impact is significant in the zone, whatever of directions, 
the wind farm's center is not entirely filled in. Table 3 presents the results of each method. 
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(a) | (b) 


Figure 7. Optimal turbine’s configuration by Grady ef al. (a) [17] current study and (b) using 30 turbines 


Table 3. Comparison results between PSO and GA 
Number of turbines Objective 10-3) _ P(kW) 
Genetic Algorithm [17] 39 1.567 17,220 
current study 39 1.376 20,187 


3.3. Case 3: Variable direction as well as speed of wind 

The third instance is the most challenging, since it is reliant on the fluctuation of both: wind speed and 
flow direction (Three speeds are taken into account in this research): 17, 12 and 8 m/s). Figure 8 shows how 
the turbines in the three-optimization methods are placed in the park boundaries to avoid the wake effectas much 
as possible. In addition, the Table 4 shows that the objective function adopted in the study is minimal in the case 
of the PSO algorithm (lower than that of GA by 0.0008403).Therefore, the output production of the wind site is 
the highest in the case of the algorithm proposed by this study (since the turbines are positioned around the 
site's perimeter, facing a free wind directly towards them).In addition, the locations of the wind turbines within 
are dispersed, It has the potential to lessen wake impacts. Table 4 shows the comparison with previously 
reported results. The suggested approach increases overall energy output while decreasing cost per unit of 
electrical generation. Figure 8 (a) and (b) depicts the optimal turbine designs. Because the most common 
occurrence of wind in the current research is between 270 and 360 degrees, as indicated in Table 1, the 
algorithm put the turbines in the upper left and lower right sides, which are linked with lesser wake impact. 


Figure 8. Optimal turbine’s configuration by Grady et al. (a) [17] current study and (b) using 39 turbines 


Table 4. Comparison results between PSO and GA 
Number of turbines _ Objective (10-3) 
Genetic Algorithm [17] 39 0.8403 
current study 39 0.776 


4. CONCLUSION 

The presented work concerned the optimization of the wind turbine’s location within the park. The goal 
of this article was to increase the energy efficiency of this wind farm while lowering the system cost. However, 
the wind device to be optimized involves several strongly coupled physical phenomena such as the wake effect. 
The study focused on the use of the swarm particle algorithm; the results are precisely improved compared to the 
same GA method that was previously used. As demonstrated throughout the study, the algorithm integrates 
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several aspects that tremendously ameliorate the performance of the previously described literature while 
simultaneously contributing to the betterment of output production with the fewest turbines. The results were 
initially compared favorably with previously published studies. The suggested algorithm’s performance and 
suitability for determining the optimal wind farm configuration have been demonstrated. To solve more complex 
cases, there is ample scope for further research in this area to investigate the problem of turbine placement 
optimization when the turbine height is different from one turbine to the other in the same sit. 
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